Self-consistent models of triaxial galaxies in MOND gravity

The Bekenstein-Milgrom gravity theory with a modified Poisson equation is tested here for the existence of triaxial equilibrium solutions. Using the nonnegative least square method, we show that self-consistent triaxial galaxies exist for baryonic models with a mild density cusp ρ ∼ Σ r . Self-consistency is achieved for a wide range of central concentrations, Σ ∼ 10 − 1000M⊙pc , representing low-to-high surface brightness galaxies. Our results demonstrate for the first time that the orbit superposition technique is fruitful for constructing galaxy models beyond Newtonian gravity, and triaxial cuspy galaxies might exist without the help of Cold dark Matter. Subject headings: galaxies: kinematics and dynamicsmethods: numerical dark matter